REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  coilection  of  information  is  estimated  to  average  1  hour  per  response,  inciuding  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  compieting  and  reviewing  this  coilection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  coliection  of  information, 
inciuding  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis 
Highway,  Suite  1 204,  Ariington,  VA  22202-4302.  Respondents  shouid  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY) 
23-03-2007 

2.  REPORT  TYPE 

Journal  Article 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

5a.  CONTRACT  NUMBER 

Effect  of  Thermal  Conductivity  on  the  Knudsen  Layer  at  Ablative  Surfaces  (Postprint) 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

L.  Pekker  (ERC,  Inc.);  M.  Keidar  (Univ.  Mich.),  J.-L.  Cambier  (AFRL/PRSA) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

50260542 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/RZSA 

10  E.  Saturn  Blvd. 

Edwards  AFB  CA  93524-7680 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-PR-ED-JA-2007- 189 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S 
ACRONYM(S) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/RZS 

5  Pollux  Drive 

Edwards  AFB  CA  93524-7048 

11.  SPONSOR/MONITOR’S 
NUMBER(S) 

AFRL-PR-ED-JA-2007- 1 89 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited  (PA  #07151  A). 

13.  SUPPLEMENTARY  NOTES 

Published  in  JOURNAL  OF  APPLIED  PHYSICS  103,  034906  2008.  ©2008  American  Institute  of  Physics 

14.  ABSTRACT 

In  this  paper  we  develop  an  analytical  model  of  Knudsen  layer  at  the  ablative  wall  taking 
into  account  the  temperature  gradient  in  the  bulk  gas.  The  region  of  validity  of  the  existing 
models  and  effect  of  the  temperature  gradient  on  the  Knudsen  layer  properties  are  calculated. 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE 
PERSON 

Dr.  Jean-Luc  Cambier 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

SAR 

7 

19b.  TELEPHONE  NUMBER 

(include  area  code) 

N/A 

standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  239.18 


Unclassified 


Unclassified 


Unclassified 


JOURNAL  OF  APPLIED  PHYSICS  103,  034906  (2008) 


Effect  of  thermal  conductivity  on  the  Knudsen  iayer  at  abiative  surfaces 

L.  Pekker,^’®^  M.  Keidar,^  and  J.-L.  Cambier^ 

^ERC  Inc.,  10  East  Saturn  Boulevard.,  Edwards  AFB,  California  93524,  USA 

^Department  of  Mechanical  and  Aerospace  Engineering,  The  George  Washington  University,  Washington, 

DC  20052,  USA 

^Air  Force  Research  Laboratory,  10  East  Saturn  Boulevard,  Edwards  AFB,  California  93524,  USA 

(Received  9  October  2007;  accepted  29  November  2007;  published  online  13  February  2008) 

In  this  article  we  develop  an  analytical  model  of  the  Knudsen  layer  at  the  ablative  wall  taking  into 
account  the  temperature  gradient  in  the  bulk  gas.  The  analysis  is  based  on  the  premise  that  the 
temperature  gradient  in  the  bulk  gas  can  be  taken  into  account  in  the  velocity  distribution  function 
at  the  outer  boundary  of  the  Knudsen  layer.  The  model  uses  a  bimodal  velocity  distribution  function 
in  the  Knudsen  layer,  which  preserves  the  laws  of  conservation  of  mass,  momentum,  and  energy  and 
converges  to  the  Chapman-Enskog  velocity  distribution  function  at  the  outer  boundary  of  the  layer. 

The  model  is  applied  to  polyethylene  ablation,  for  which  two  cases  are  considered:  (a)  the  ablation 
process  is  due  to  pure  heat  conduction  to  the  surface,  with  no  external  heating  of  the  ablated  surface, 
and  (b)  the  ablation  is  due  to  both  the  thermal  conduction  and  an  external  heating  of  the  surface,  e.g., 
vaporization  of  a  metal  exposed  to  laser  radiation.  The  region  of  validity  of  the  existing  models  and 
effect  of  the  temperature  gradient  on  the  Knudsen  layer  properties  are  calculated.  ©  2008  American 
Institute  of  Physics.  [DOI:  10.1063/1.2838210] 


I.  INTRODUCTION 

The  physics  of  the  Knudsen  layer,  the  nonequilibrium 
(kinetic)  layer,  formed  near  the  vaporizing  (ablative)  surface 
is  of  great  interest  for  a  number  of  applications,  such  as 
capillary  discharges, plasma  thrusters,^’"^  high-pressure 
discharges,^  vacuum  arcs,^  electroguns,^  and  laser  ablation.* 

Anisimov^  was  the  first  to  consider  details  of  the  vapor¬ 
ization  process  for  a  case  of  vaporization  of  a  metal  exposed 
to  laser  radiation.  He  used  a  bimodal  velocity  distribution 
function  in  the  kinetic  layer,  assuming  no  absorption  of  laser 
radiation  in  the  ablated  gas.  The  primary  result  of  his  work 
was  the  calculation  of  the  maximal  flux  of  returned  atoms  to 
the  evaporating  surface,  which  was  found  to  be  about  18%  of 
the  flux  of  vaporized  atoms.  This  result  was  obtained  under 
the  assumption  that  the  atom  flow  velocity  is  equal  to  the 
sound  velocity  at  the  external  boundary  of  the  Knudsen  layer 
and  the  temperature  of  the  gas  in  the  equilibrium  region  (be¬ 
yond  the  Knudsen  layer)  is  constant,  i.e.,  no  conductive  heat 
flux  to  the  ablative  wall  surface. 

However,  in  many  physical  situations,  the  vapor  leaving 
the  nonequilibrium  layer  cannot  be  described  by  using  a 
speed  of  sound  approximation.  For  example,  in  ablative  cap¬ 
illary  discharges,  the  gas  motion  in  the  capillary  chamber  is 
not  “free”;  it  is  restricted  by  the  capillary  wall,  leading  to  a 
more  dense  gas  (plasma)  in  the  discharge  volume  and  there¬ 
fore,  larger  back  flux  to  evaporating  surface  and  smaller  flow 
velocity  at  the  outer  boundary  of  the  Knudsen  layer. 
Ytrehus*®  has  used  the  Anisimov  and  Ansatz  bimodal  veloc¬ 
ity  distribution  functions  in  the  Knudsen  layer  to  study  the 
effect  of  bulk  gas  pressure  on  downstream  vapor  flow  (half¬ 
space  evaporation  problem).  He  has  calculated  the  density 
and  temperature  jumps  over  the  Knudsen  layer,  the  evapora- 


“^Electronic  mail:  leonid.pekker.ctr@edwards.af.mil. 
0021  -8979/2008/1 03(3)/034906/6/$23.00 


tion  mass  flux  and  other  parameters  of  the  Knudsen  layer  as 
functions  of  the  ratio  of  the  equilibrium  vapor  pressure  to  the 
gas  pressure  at  the  outer  boundary  of  the  Knudsen  layer  (gas 
bulk  pressure).  In  addition,  Ytrehus  has  demonstrated  that  his 
analytical  results  are  in  substantial  agreement  with  the  ex¬ 
perimental  finding,  direct  simulation  Monte  Carlo  (DMCS) 
and  numerical  solutions  of  the  Boltzmann  equations.  He  has 
also  shown  that  the  differences  between  the  analytical  solu¬ 
tions  using  the  Anisimov  approximation  and  the  Ansatz 
(more  sophisticated)  velocity  distribution  function  are  very 
small. 

Beilis  ’  was  the  first  to  consider  ablation  into  a  dense 
plasma.  He  studied  the  case  of  metal  vaporization  into  dis¬ 
charge  plasmas  in  a  vacuum  arc  cathode  spot.  He  concluded 
that  the  parameters  at  the  outer  boundary  of  the  Knudsen 
layer  are  close  to  their  equilibrium  values  and  that  the  veloc¬ 
ity  at  the  outer  boundary  of  the  kinetic  layer  is  much  smaller 
than  the  sound  velocity. 

Later  these  models  were  applied  for  the  case  of  dielectric 
ablation  into  the  discharge  plasma  in  the  capillary  discharge 
conditions**’'"^  and  for  the  case  of  strong  plasma 
acceleration.**  All  those  analytical  models  neglected  the  con¬ 
ductive  heat  flux  to  the  ablative  surface.  This  can  be  signifi¬ 
cant  because  the  temperature  in  the  plasma  core  is  assumed 
in  the  models  to  be  much  greater  than  the  temperature  of  the 
ablative  surface.  In  particular,  neglecting  the  conductive  heat 
flux  results  in  the  calculated  gas  temperature  at  the  outer 
boundary  of  the  Knudsen  layer  to  appear  to  be  smaller  than 
the  temperature  of  the  evaporating  surface.  This  conse¬ 
quently  leads  to  the  heat  flux  through  the  Knudsen  Layer 
being  directed  upward  to  the  plasma  chamber.  Therefore,  ne¬ 
glecting  the  conductive  heat  transfer  in  the  Knudsen  layer 
leads  to  an  inconsistency  in  all  the  models  where  the  gas 
(plasma)  temperature  is  larger  than  the  surface  temperature 
of  the  ablative  wall. 
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However,  it  is  worth  noting,  that  in  the  case  where  the 
external  heat  flux  to  the  ablating  surface  is  larger  than  the 
conduction  heat  flux,  the  temperature  of  the  ablating  surface 
might  be  larger  than  the  gas  temperature  at  the  outer  bound¬ 
ary  of  the  Knudsen  layer.  This  is  the  case  for  the  example  in 
Ref.  9  where  an  externally  applied  laser  radiation  source 
heats  the  ablating  surface,  but  for  which  the  gas  (plasma)  is 
transparent. 

Recently  many  numerical  models  of  evaporation  pro¬ 
cesses  were  developed  based  on  Monte  Carlo  simulation*^ 
and  numerical  solutions  of  the  Boltzmann  equations,*^  de¬ 
scribing  the  kinetic  layer  without  any  a  prior  approximation 
of  gas  velocity  function  distribution  in  that  layer.  However, 
as  Rose  has  demonstrated  the  models  using  his  and  Anisi¬ 
mov’s  bimodal  velocity  distribution  functions  in  the  kinetic 
layer  give  results  virtually  coincident  with  numerical  solu¬ 
tions  of  the  Boltzmann-Krook-Welander  equation  for  evapo¬ 
ration  of  a  monatomic  substance  with  condensation  coeffi¬ 
cient  equal  to  unity.  As  shown  in  Refs.  10  and  16-18 
analytical  models  also  are  in  good  agreement  with  Monte 
Carlo  simulation. 

Ideally  Monte  Carlo  simulations  should  be  able  to  self- 
consistently  describe  the  conductive  heat  flux  to  the  ablating 
surface.  However,  this  will  require  extending  the  analysis 
beyond  the  Knudsen  layer  region,  making  it  computationally 
intensive.  Thus,  improving  the  analytical  models  by  includ¬ 
ing  heat  conduction  into  consideration  is  an  important  step  in 
developing  practical  (computationally  efficient)  solutions  for 
modeling  of  evaporation  processes  and  plasma  discharges 
coupled  to  ablative  processes,  and  in  improving  our  physical 
understanding  of  the  Knudsen  layer. 

We  would  also  like  to  point  out  the  recent  article  by 
Bond  and  Struchtrup,  in  which  the  authors  have  included 
the  conduction  heat  flux  in  their  analytical  model  of  water 
evaporation.  This  is  a  generalized  Hertz-Knudsen  model  of 
the  Knudsen  layer,  which  uses  the  condensation  and  accom¬ 
modation  coefficients  (probabilities)  at  the  evaporating  sur¬ 
face  to  take  into  account  back-flux  effect  in  the  Knudsen 
layer.  However,  these  coefficients  can  have  a  complex  de¬ 
pendence  on  the  evaporation  conditions  such  as  vapor  pres¬ 
sure,  temperature,  surface  conditions,  incidence  angle,  etc., 
and  are  usually  determined  experimentally  or  by  fitting  the 
model  with  experimental  data.  In  the  case  of  dielectric  abla¬ 
tion  into  plasma,  where  the  ablation  parameters  vary  among 
a  wide  range  of  temperatures  and  pressures,  these  coeffi¬ 
cients  are  not  experimentally  measured.  It  is  worth  noting 
that  the  classical  Hertz-Knudsen  model  and  all  its  genera¬ 
tions  assumes  no  collisions  in  the  Knudsen  layer  (i.e.,  does 
not  satisfy  conservation  of  momentum  through  the  Knudsen 
layer)  preventing  any  “relaxation”  in  the  kinetic  (Knudsen) 
layer,  where  the  velocity  distribution  function  at  the  ablative 
surface  has  to  relax  (converge)  to  the  bulk  gas  distribution 
function  at  the  outer  boundary  of  the  Knudsen  layer.  The 
analytical  models  using  bimodal  velocity  distribution  func¬ 
tions,  such  as  the  one  described  here,  take  into  account  the 
collisions  in  the  Knudsen  layer,  satisfy  the  conservation  of 
momentum  through  the  Knudsen  layer,  self-consistently  cal¬ 
culate  the  back  flux  and  implicitly  assume  a  condensation 
coefficient  of  unity.  As  has  been  mentioned  earlier,  they  de¬ 


scribe  the  ablation  process  with  reasonable  approximation, 
and  the  effect  of  variable  condensation  coefficient  is  not  con¬ 
sidered  there. 

In  this  article  we  develop  an  analytical  model  of  the 
Knudsen  layer  by  considering  an  appropriate  boundary  con¬ 
dition  in  the  kinetic  formulation  that  takes  into  account  the 
gas  temperature  gradient  at  a  flat  gas-wall  interface.  The 
region  of  validity  of  the  existing  models  and  the  effect  of  the 
temperature  gradient  on  the  Knudsen  layer  properties  are  cal¬ 
culated.  The  main  impetus  of  this  article  is  to  study  the  effect 
of  the  thermal  conductivity  on  the  Knudsen  layer  formed 
near  the  ablated  surface.  This  analysis  is  based  on  the 
premise  that  thermal  conductivity  (the  temperature  gradient) 
in  the  gas  bulk  can  be  taken  into  account  in  the  velocity 
distribution  function  at  the  outer  boundary  of  the  Knudsen 
layer.  In  this  article  we  use  such  a  function  obtained  by 
Chapman-Enskog  expansion  method  for  solving  the  Boltz- 
mann  equation  based  on  the  assumption  that  the  molecular 
mean-free  path  is  much  smaller  than  characteristic  scale  of 
the  temperature  change.  Thus,  our  model  is  limited  to  rela¬ 
tively  small  values  of  the  temperature  gradients. 


II.  MODEL  FORMULATION 

Following  Anisimov’s  method,®  let  us  write  the  velocity 
distribution  function  in  the  kinetic  layer  with  the  evaporating 
surface  in  the  following  form.  Fig.  1 

/(x,V)  =  Six)  •/,(¥)  +  [1  -  Six)]  •/„(¥),  (1) 
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Here  fi,  is  the  velocity  function  distribution  at  the  inner 
boundary  of  Knudsen  layer  (at  the  ablative  surface)  with 
Maxwellian  vaporization  function  for  y^>0  and  a  shifted 
“back  flux”  Maxwellian  function  of  the  particle  for  V^<0 
describing  the  particles  incoming  to  the  surface  from  the  gas, 
where  the  x  axis  is  normal  directed  to  the  wall  from  the  gas 
chamber;  /„  is  the  Chapman-Fnskog  velocity  distribution 
function  at  the  outer  boundary  of  the  Knudsen  layer  that 
takes  into  account  the  temperature  gradient  and  directed 
velocity  above  the  Knudsen  layer,  as  shown  in  Fig.  1; 
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ni,  Ti,  f 


Gas  bulk»X 


”^0' 


Knudsen  layer~A,|^fp 


FIG.  1.  Schematic  representation  of  the  layer  structure  near  the  ablative 
surface. 


kTQ=mVj/2  is  the  temperature  of  the  ablated  surface;  v  is 
the  collision  frequency  depending  on  the  temperature  and 
density  of  the  gas;  is  an  unknown  function  that  satished 
the  conditions  <5(0)=1  and  ^oo)  =  0.  The  parameter  /3  is  an 
unknown  variable  that  must  be  obtained  in  the  solution;  it 
represents  essentially  nonequilibrium  effects  caused  by  col¬ 
lisions  in  the  Knudsen  layer.  The  number  densities  are  nor¬ 
malized  on  the  equilibrium  vapor  number  density  Hq  corre¬ 
sponding  to  the  surface  temperature  Tq  and  all  velocities  are 
normalized  on  Vj-.  The  expression  for  takes  into  account 
the  fact  that  the  vaporized  atoms  have  a  Maxwellian  distri- 

23 

bution  at  a  temperature  equal  to  the  surface  temperature 
and  also  assumes  that  the  number  density  of  the  evaporated 
atoms  is  equal  to  half  of  the  equilibrium  vapor  number  den¬ 
sity  at  this  surface  temperature,  i.e.,  reasonable  approxima¬ 
tions  used  in  all  previous  bimodal  velocity  distribution  func¬ 
tions  models.'^  It  should  be  noted  that  Anisimov’s  model^ 
(as  well  as  all  other  existing  models)  employed  shifted  Max¬ 
wellian  function  at  the  outer  boundary  of  Knudsen  layer  ig¬ 
noring,  as  mentioned  previously,  the  conduction  heat  transfer 
to  the  ablative  surface,  i.e.,  temperature  gradient  at  the  edge 
of  the  Knudsen  layer. 

Assuming  the  conservation  laws  of  mass,  momentum, 
and  energy  hold  at  all  times  within  the  discontinuity  region, 
through  the  Knudsen  layer,  as  it  has  been  assumed  in  all 
previous  models^”^^  (quasisteady  state  approximation  within 
the  Knudsen  layer),  the  following  integrals  are  defined: 

^+00  rco  r+00 

Cl  =  I  dVA  dVyl  fV,^dV^  =  n^u,  (5) 


M^  =  mnoVjCi, 

(5') 

Px  =  mnoVlC2, 

(6') 

vl 

Ex  =  mnQ—C2„ 

(7') 

where  consists  of  the  two  parts:  the  conduction  heat  flux 
and  the  enthalpy  flux  of  the  gas  moving  with  a  directed  ve¬ 
locity,  Eq.  (7).  Taking  into  account  that  integrals  Ci,  C2,  and 
C3  are  preserved  through  the  Knudsen  layer  and  they  should 
be  independent  on  5(x),  we  obtain  the  following  equations 
corresponding  to  Ci,  C2,  and  C3: 


2^  =  C,-n,(3 
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Equations  (5),  (5'),  (6),  (6'),  (8),  and  (9)  are  identical  to  the 
corresponding  mass  and  momentum  conservation  equations 
obtained  in  Refs.  11-18,  whereas  Eqs.  (7),  (7'),  and  (10) 
differ  from  the  corresponding  energy  conservation 
equations  by  the  temperature  gradient  term,  which  is  re¬ 
sponsible  for  conduction  heat  flux  to  the  ablative  surface. 

Let  us  introduce  a  thermal  conduction  parameter  Tj 


VjVi  <7(ln  T) 
V  dx 


(11) 


where  knifp=(  Vj-  Vi)/  v  is  the  gas  mean-free-path  at  the  outer 
boundary  of  the  kinetic  layer  and  Sxj=\^d{\n  T)l is  the 
characteristic  gradient  length.  Condition  (11)  is  needed  for 
the  Chapman-Enskog  expansion  method  and  Eq.  (7)  to  be 
valid,  as  explained  earlier.  In  the  case  of  small  u,  Eqs. 
(8)-(10)  can  be  simplihed  to  the  following  form: 
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where  values  of  Ci,  C2,  and  C3  obtained  at  the  outer  bound¬ 
ary  of  the  Knudsen  layer,  where  S  is  equal  to  zero,  and  the 
mass,  momentum,  and  energy  fluxes  are 


As  we  can  see,  for  u  =  0  (in  which  case  there  is  no  abla¬ 
tion)  and  for  r7->0,  the  temperature  and  density  at  the  outer 
edge  of  kinetic  layer  are  correspondingly  larger  and  smaller 
than  1 ;  it  is  as  expected,  as  the  gas  bulk  region  has  a  higher 
temperature  than  the  wall  surface  and  the  total  conduction 
heat  flux  is  assumed  fo  be  directed  to  the  wall  (t7->0).  These 
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FIG.  2.  Normalized  u,  Ti  =  Vl,  and  back  flux  at  the  ablative  surface  vs  Tj-. 


FIG.  3.  Comparison  of  exact  and  approximate  solutions.  Polyethylene  with 
a  surface  temperature  of  800  K. 


temperature  jumps  between  the  solid  surface  and  the  Max¬ 
wellian  gas  at  the  outer  edge  of  the  kinetic  layer  is  a  very 
well  known  phenomenon  in  gas  dynamics,  see,  e.g..  Ref.  24. 

Assuming  that  the  ablation  process  is  due  to  only  the  gas 
thermal  conduction  (no  external  heat  source  is  applied  to  the 
ablative  surface),  the  boundary  condition  at  the  ablative  sur¬ 
face  can  be  written  as 

-£x-'^wall-^=M,-d>,,p,  (15) 

where  is  thermal  conductivity  of  the  wall,  is 

evaporation  heat  of  the  wall  material,  -E^  is  the  total  energy 
flux  through  the  Knudsen  layer  incoming  into  the  ablative 
wall,  Eq.  (7'),  and  is  the  mass  ablation  rate,  Eq.  (5');  the 
negative  sign  in  front  of  is  due  to  the  x  axis  being  directed 
from  the  wall  into  the  gas  chamber.  Assuming  no  heat  loss  in 
the  bulk  of  the  wall,  the  boundary  condition  at  the  ablative 
surface  is  reduced  to: 

-£,  =  M,-<Dvap.  (16) 

This  can  be  the  case  of  a  “boiling  wall,”  where  all  heat 
incoming  into  the  wall  is  spent  on  vaporization.  Thus,  for  a 
given  ablative  surface  temperature  and  corresponding  equi¬ 
librium  vapor  pressure  and  number  density,  and  given  heat 
conduction  parameter  Tj,  Eqs.  (8)-(10)  and  (15)  can  be 
solved  relative  to  variables  nj,  Vi,  u,  and  13.  The  total  abla¬ 
tion  rate  and  the  heat  flux  to  the  ablative  wall  [Eqs.  (5')  and 
(7')]  can  then  be  computed.  The  dependences  of  rii,  Vi,  and 
back  flux  (the  total  flux  of  particles  incoming  to  the  ablative 
surface  from  the  gas)  on  the  Tj  and  for  a  specific  example  of 
polyethylene  wall  are  presented  in  the  following. 

III.  RESULTS 

Figure  2  shows  the  calculated  parameters  of  the  Knudsen 
layer  versus  Tj  for  the  case  of  thermal  conduction  heating  of 
the  ablative  polyethylene  wall  with  the  polyethylene  surface 
temperatures  of  650  and  800  K,  Eqs.  (8)-(10)  and  (16);  the 
evaporation  heat  has  been  taken  as  3.6X10®  (J/kg)  and 
equilibrium  vapor  pressure  equal  as  P=  10®exp[5565.22(l/ 
453-1/7’)],^®  where  the  pressure  is  in  pascal  and  tempera¬ 


ture  is  in  kelvin.  One  can  see  that  u,  T,  back  flux,  and  are 
weak  functions  of  temperatures  in  this  temperature  region, 
although  their  equilibrium  vapor  pressures  differ  almost  by 
five  times.  As  it  has  been  mentioned  earlier,  we  cannot  ex¬ 
tend  the  obtained  results  of  Fig.  2  for  higher  Tj,  because  the 
Chapman-Enskog  expansion  method  is  valid  only  for  Tj-<^  1 . 
In  this  limit,  Eq.  (16)  can  be  simplified  to  the  following 
form: 


-  u 


kTo  I 


(17) 


where  m  is  the  average  atomic  mass  of  polyethylene  compo¬ 
sition  and  Tq  is  the  polyethylene  surface  temperature.  Ex¬ 
pressing  M  as  a  function  of  Tj,  Eq.  (17),  and  substituting  it 
into  Eqs.  (12)  and  (14)  and  then  into  the  following  equation 
for  normalized  back  flux  at  the  ablative  wall: 

I  -  unil'fn— 1  —  u2'[tt,  (18) 

yields  the  explicit  relationships  between  «[,  T=v\,  u, 
and  Tj.  Comparison  of  this  approximate  solution  and  the 
“exact”  solution  obtained  from  the  solution  of  Eqs.  (8)-(10) 
is  shown  in  Fig.  3,  leading  to  the  conclusion  that  approxi¬ 
mate  solution  gives  satisfactory  results  for  Tj.<0.05. 

We  also  would  like  to  point  out  that  in  the  case  of  ther¬ 
mal  conduction,  u  is  small  and  cannot  reach  sonic  condition. 
Otherwise,  E^  would  be  positive,  Eqs.  (7)  and  (7'),  meaning 
that  the  total  energy  flux  would  be  directed  not  to  the  abla¬ 
tive  surface  but  upward  into  the  gas  chamber,  contradicting 
the  model  assumptions. 

In  the  case  where  the  ablating  surface  is  heated  by  an 
additional  (external)  heat  source,  e.g.,  by  laser  radiation,^  Eq. 
(16)  can  be  rewritten  as 

-£x-£ext=A7,‘I^vap,  (19) 

where  is  an  external  heat  flux  to  the  ablative  surface.  It  is 
worth  noting  that  in  the  case  where  thermal  conduction  in  the 
wall  is  not  equal  to  zero,  see  Eq.  (15),  the  E^.^i  is  the  net 
external  heat  flux  at  the  wall  equal  to  the  external  heat  flux  to 
the  wall  “above”  the  surface  minus  conduction  heat  into  the 
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FIG.  4.  Normalized  thermal  and  directed  velocities  at  the  outer  boundary  of 
the  Knudsen  layer  as  functions  of  ratio  of  external  heat  flux  to  the  total  heat 
flux  at  the  ablating  polyethylene  surface. 

wall  “behind”  the  surface.  With  an  increase  in  the  abla¬ 
tion  rate  increases  and,  if  becomes  much  larger  than 
conduction  heat  flux,  we  may  drop  the  conduction  heat  term 
from  the  velocity  distribution  function  [Eq.  (3)],  recovering 
the  previous  models.®  Figure  4  shows  Vi  and  u  as  func¬ 
tions  of  the  ratio  of  an  external  heat  flux  to  the  total  heat  flux, 
g=£ext/(-E^ext+£'thermai)  fot  the  ablative  polyethylene  wall 
with  the  polyethylene  surface  temperatures  of  800  K  and  for 
r7'=0.01,  0.1,  and  0.3;  the  calculations  have  been  performed 
up  to  sonic  conditions  with  7=5/3.  As  one  can  see  with  a 
decrease  in  and  an  increase  in  q,  the  distributions  of  Vi 
and  u  are  converging  to  the  case  of  t=0. 

The  effect  of  heat  conduction  in  the  case  of  an  external 
heat  flux  can  also  be  calculated,  following,  as  functions 
of  a=ulVi  for  different  Tj,  Eqs.  (8)-(10);  here,  the  variation 
of  a  indicates  the  magnitude  of  the  external  heat  flux  com¬ 
pared  to  the  conduction  heat  flux  (as  a  and  the  flow  velocity 
increase  up  to  sonic  condition,  there  can  be  no  heat  conduc¬ 
tion  to  the  wall).  The  results  of  these  calculations  are  shown 
in  Fig.  5.  As  one  can  see  at  small  a  the  effect  of  thermal 


FIG.  5.  Normalized  temperature  and  density  at  the  outer  boundary  of  Knud¬ 
sen  layer  as  a  function  of  a=u/  Vi  with  temperature  gradient  at  the  Knudsen 
layer  edge  as  a  parameter. 


conduction  is  important,  and  with  an  increase  in  a,  the  effect 
of  thermal  conduction  decrease,  increases  and  changes 
sign  [becomes  positive  meaning  that  the  total  energy  flux  is 
directed  upward,  Fig.  1,  into  the  gas  chamber,  see  Eqs.  (7) 
and  (7')]  and  the  temperature  and  density  plots  are  converg¬ 
ing  to  the  case  of  T7-=0.  As  one  can  see,  Ti  and  rii  decreases 
with  an  increase  in  a  (with  an  increase  in  directed  velocity 
m).  Fig.  5,  that  can  be  explained  by  expansion  of  the  dense 
ablated  gas  stream  in  less  dense  gas  surrounding.  It  is  worth 
noting,  that  such  a  decrease  in  T^  and  Wj  leads  to  a  decrease 
in  the  thermal  heat  conduction  at  a  given  Tj,  see  Eqs.  (7)  and 
(11). 


IV.  CONCLUDING  REMARKS 

In  summary,  a  model  of  the  Knudsen  layer  near  the  ab¬ 
lated  surface  in  the  case  of  surface  heating  by  the  adjacent 
gas,  by  virtue  of  thermal  conductivity,  was  developed.  Pre¬ 
viously  existing  models  were  not  able  to  describe  this  physi¬ 
cal  situation  due  to  neglecting  effect  of  thermal  conductivity 
on  the  velocity  distribution  function  in  the  gas  and  thus  lead¬ 
ing  to  temperature  gradient  directed  outwards  the  surface.  In 
contrast,  the  developed  model  predicted  existence  of  the  tem¬ 
perature  gradient  toward  the  surface.  It  should  be  pointed  out 
that  this  model  is  limited  to  relatively  small  temperature  gra¬ 
dients  due  to  the  limitation  of  the  Chapman-Enskog  expan¬ 
sion  for  solving  the  Boltzmann  equation.  In  the  case  of  a 
larger  temperature  gradient  a  more  rigorous  model  is  re¬ 
quired  and  only  numerical  simulations  such  as  DSMC  would 
be  able  to  solve  the  problem.  Finally,  the  model  can  be  also 
verified  experimentally,  e.g.,  in  liquid-vapor  experiments, 
where  the  surface  temperature,  the  mass  flux,  and  the  bulk 
gas  pressure  are  measured  as  in  Ref.  26.  Unfortunately,  we 
cannot  use  this  work,  because  their  experiments  have  been 
performed  on  “bending”  liquid-vapor  interface.  As  shown  in 
this  article  and  in  Ref.  20,  the  measured  high  temperature 
jump  in  this  experiment  is,  probably,  due  to  very  small  cur¬ 
vature  radius  of  the  vapor-liquid  interface.  To  the  best  of  our 
knowledge  there  currently  is  no  experimental  data  that  can 
allow  direct  experimental  verification  of  the  model,  and  com¬ 
parison  with  direct  numerical  simulations  will  play  a  similar 
role  until  a  more  complete  physical  model  can  be  developed. 
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